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By applying the renormalization group method to two-coupled chains in the Tomonaga 
model, the role of interchain hopping has been studied in the entire energy region. The energy 
for a crossover from the perturbational regime to the relevant regime becomes smaller than 
that of the interchain hopping due to one-dimensional fluctuations of the mutual interaction. 
From the calculation of response functions for charge density waves and superconducting 
states, the phase diagram of dominant and subdominant states has been obtained in the 
plane of mutual interactions with fixed energy. 

§1. Introduction 

Organic conductors, & which exhibit low dimensional fluctuations, have been 
studied by using the model of quasi-one-dimensional systems, where the role of the 
interchain hopping plays an important role. 

As a basic model for such a system, two-coupled chains of interacting electrons 
has been explored. Even in the simple case of spinless fermions described by the 
Tomonaga model, there exist significant properties involved in connecting a one- 
dimensional system with quasi-one dimensional system. Ef _ Q) The ground state has 
been calculated in the limit of low energy. The relevant behavior of the interchain 
hopping leads to a state which differs essentially from that in the absence of hopping. 
Although the limiting state for the two-coupled chains of the Tomonaga model is well 
known, it is not yet clear how the state at finite temperature (or energy) is determined 
in the presence of both the interchain hopping and the mutual interactions. 

The phase diagram in the plane of the intrachain interaction (#2) and the inter- 
chain interaction (g' 2 ) is shown in Fig. 1(a) (Fig. 1(b)). 0) This represents the result 
in the absence (presence) of the interchain hopping. In Fig. 1 (a) , there are two kinds 
of dominant states with in-phase pairing and out-of-phase pairing which are degener- 
ate due to the absence of interchain hoping. Such a state corresponds to the state in 
the high energy limit, where the effect of the interchain hopping can be disregarded. 
The state in Fig. 1(b) is obtained in the limit of zero energy (or temperature). The 
interchain hopping plays a crucial role in the sense that the degeneracy for in-phase 
and out-of-phase pairings is removed and the other pairing state becomes subdomi- 
nant due to the relevance of the interchain hopping. The subdominant state which 
is shown in the parenthesis could become meaningful for more complicated models, 
e.g., a model with a Fermi surface which exhibits an incomplete nesting condition 
in the case | £72 1 > |$2 1- Thus, it is of interest to examine the state with finite energy, 
which is expected to exhibit a phase diagram between Figs. 1(a) and (b). 
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(a) (b) 

Fig. 1. [Phase diagram in the plane of g 2 = 72/(27t«p) and g' 2 = r y' 2 /(2TTVF) for t = (a) and t ^ 
(b),EP where C and 5* denote the CDW and SC state, respectively, and || (_L) and out (in) denote 
the paring state for the intrachain (interchain) and out of phase (in phase) ordering. For t = 
(a), there are two kinds of states as the dominant states, while for t 7^ (b), the subdominant 
state is distinguished, as is shown in the parenthesis. 

In deriving Fig. 1(b), it has been assumed that one-dimensional behavior van- 
ishes and the relevant behavior of the interchain hopping begins when the energy (or 
temperature) becomes smaller than the hopping energy. In a quasi-one dimensional 
system, Bourbonnais u has maintained that one-dimensional fluctuations still exist 
at temperatures lower than the hopping energy and that the temperature for the 
crossover becomes smaller than the hopping energy due to the suppression of the 
density of state around the Fermi energy. u> There is no explicit calculation which 
shows the reduction of the crossover temperature, even for two-coupled chains. 

In the present paper, we study these problems in two-coupled chains of the 
Tomonaga model by applying the renormalization group method to the bosonized 
Hamiltonian. In §2, the Hamiltonian and possible order parameters are represented 
in terms of the phase variables. The scaling equations for coupling constants and 
response functions are given in §3. In §4, a phase diagram with a fixed energy is 
examined, and the crossover energy is evaluated as a function of mutual interactions. 
Section 5 is devoted to summary and discussion. 

§2. Model and order parameter 

We consider two-coupled chains of the Tomonaga model, where a one-dimensional 
chain consists of spinless fermions with intrachain forward scattering, interchain for- 
ward scattering, and interchain hopping. The Hamiltonian with a chain of length L 
is given as 



H 



E 

p.t 



dx V>p j(aO v F (~ipd x - k F ) i/> P)i {x) 
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-*E / dx{-ipl- L (x)'ip Pj 2(x) +h.c.j 
v J 

+ yE/ dx ^l,i( x )^- P ,i( x )^-vA x )^pA x ) 

p,i 

+ y E J dx {44 x ^-pA x ^-p^ x )M x ) + h - c -} > ( 24 ) 

where % is the fermion creation operator. The sign p = + (— ) expresses the right 
(left) moving state, and i = 1 (2) denotes the indices of the chain, where vp and kp 
are the Fermi velocity and the Fermi wave vector, respectively. We consider the case 
in which the energy of interchain hopping, t, is much smaller than the Fermi energy, 
vpkp, and the matrix element of the intrachain (interchain) forward scattering, 72 
(72), is smaller than 2nvp. 

In Eq. (2T) the first and second terms can be diagonalized by use of the trans- 
formation 

1>pA x ) = {-^pAx) + V P ,2(s))/\/2. (2-2) 
The index a = ± denotes that for the diagonalized band, where the new Fermi wave 
vector is given by kp a = kp — at/vp for the cr-band. By applying the bosoniza- 
tion method to the respective band, Eq. (2-1) is rewritten in terms of the phase 
Hamiltonian as H = Tig + ?^,Ef'0) where 

n e = ^J dx^(d x e + f + K 6 (d x e.) 2 Y (2-3) 



n * = ti J dx ) JT {dx<p+) + K *^-) 



+ ^^2 / dx {^20+ cos (2(/> + - -^-x\ + g 2( p- cos (2<£_) j . (2-4) 

Equations (2-3) and (2-4) express the Hamiltonian of the total fluctuation and that 
of the transverse fluctuation, respectively, where u = vp J 1 — (g 2 + g'2) 2 , Kg = 



{1 - (92 + 3 2 )}/{ 1 + (92 +92)}, K 4> = 1, 92<t>+ = -(32-52), 320- =92-92 
and 72 = 2irvpg2 (72 = ^.irvpg^). The quantity a , which is of the order of kp, is 
the cutoff for large wave vectors. The phase variables 6± and <j>± in Eqs. (2-3) and 
(2-4) are defined as 

"±0»0 = ~ E ^e~-M^ £ (4,+,^k +q , + ,a ± i/H,-,) , (2-5) 

TT'l 

q^O qL k,a 

where afc iPiCr = (1/a/Z) / dx e~ ik ^jp Pta (x). By use of Eqs. (2-5) and (2-6), the field 
operator ip Pta (x) is expressed asE^ 

^pA x ) = at — ex P [ipkpaX + pO P)(J + iirEpa] = ip' (x)exp[iir~ p o -], (2-7) 
\/27ra 



M*) = "E ^ e " f M+ ^E^ (4,+,^+*,+,* ±a{_ ta a k+q ,^) , (2-6) 
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where & Pja = i[6 + + pO^ + a(<p + + p<p-)]/2. The phase factor ifS P]ff , which is 
introduced for the fermion operator, is given by = 0, = —-,+ = 

N+,+ + N+- and ~_ _ = iV +i+ + + where TV^ = / dx V P , ff (x) 

is the number operator for the fermion. 

We examine the state with order parameters for the charge density wave (CDW) 
state and the superconducting (SC) state. In terms of phase variables, order param- 
eters corresponding to Fig. 1(a) are expressed as follows. 
(I) Order parameters for CDW with intrachain and out of (in) phase pairing are 
given by 



1 1, out 



c 



1 1, in 



<j 



— i2pk^x 

e -ip8+ 

7TQ 



-z sm 0_ 
2t 

COS X 



(2-8) 



(II) Order parameters for CDW with interchain and out of (in) phase pairing are 
given by 



C_L,out 
C_|_,in 



( Wp^-P? - Vp,2^-P,l 
\ ^,1^,2+^,2^,1 

/ E<,^-p,-a \ 



(7 



J2 a ^p,^-p,a 



J 



-i2pkpx 



ira 



COS (j)- 



-ipsm 



2t 

Vf 



(2-9) 



(III) Order parameters for the superconducting state with intrachain and in (out of) 
phase pairing are given by 



% out 



1pp,l1p-p,l + 1pp,2lp-p,2 

(7 



Y^4>p,ai>-p,a \ 

(7 



Crossover in Two- Coupled Chains of Tomonaga Model 



5 



/ COS0- \ 

= — e~ tpe ~ . . / , 2t \ . (2-10) 
7ra y zpsin I #> + — — x I y 

(IV) Order parameters for the superconducting state with interchain and in (out of) 
phase pairing are given by 

/ / / , / / \ / -J2 a ^p,^-p,° 




^V^p^-p^ 



1 



z sm ( 



= — e ' ipd ~ (a 21 \ • ( 2 ' n ) 

7ra y cos I + - —a; I y 

The symbol || (_L) represents the pairing in the chain (between two chains), and the 
symbol in (out) denotes pairing with the same (opposite) sign. The right arrow in 
the second line of Eqs. (2-8)~(2-ll) denotes the process obtained from Eq. ( |2-7p , 
where a numerical factor with absolute value of unity is discarded. In Eq. (2-4) and 
Eqs. (2-8)~(2Tl), the quantity N+ + — iV_ + is taken as an even integer, where N PttT 
is an eigenvalue of N Pi(T . We note that the symmetry of the interactions is associated 
with the four kinds of order parameters of Eqs. (2-8)~(2Tl). Since a quarter of 
the 52-52 pl ane i n Fig. 1(a) is sufficient to understand all the states, we precisely 
examine the states in the region 52 > Ifi^l- 

§3. Renormalization group method 

Several kinds of order parameters introduced in the previous section are exam- 
ined by calculating the response functions given by 

R a (x,t) = (t t A^(x,t)A(0,0)\ ■ 2(vra) 2 

= R a (x,t) ■ R* a (x,t), (3-1) 

where A is the operator given by Eqs. (2-8)~(2Tl). From Eq. (2-3), R a (x,t) is 
calculated straightforwardly as 

/ \ K ^ 

R a (x,t) = ( Tre ^) e -^o,o)\ = » , ( 3 .2) 



where the sign + (— ) corresponds to the CDW (SC) state. Note that Eq. 
depends only on r(= x 1 + (vpr) 2 ). From Eqs. (2-8)~(2Tl), the quantity R a (x,t) 
corresponding to the 0-field is expressed as 

R* a (x,t) = R^(x,t) ■ f A (x), (3-3) 
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where 



( f? 

- rl sin< 



B* A {Z,T) 



R. 



'COS 0- 



R 



sin 0+ — 



R, 



'COS <f>+ — 



2 (T T sin 0_ (x, r) sin 0_ (0, 0)) 

for A = C|| iOUt ,5_L )in , 
2 (T T cos 0_ (x, r) cos 0_ (0, 0)) 

for A = Cj_ iOUt ,5j| )in , 
2 (T T sin (0+(x, r) + 2tx/v F ) sin 0+(O, 0)) 

cos (2tx/vp) 

for ^ = Cj_ iin , Sp )OUt , 
2 (T r cos (<p+(x, t) + 2tx/u F ) cos </>+(0, 0)) 

cos {2tx/v-p) 

for A = Cm in , Si out 



(3-4) 



and the factor /a(x) is defined as Ja{x) = e l2pk¥X for A = Ci out and Cj_ )OU t, Ja{x) = 
e i2 P k F x cog (^x/vf) for A = C|| jin and C ± , in , / A (x) = 1 for A = S|| )in and S ± , in , and 
/a(x) = cos (2tx/v F ) for A = S|| )OUt and S_L )OUt . 

By using the renormalization group method with the assumption of invariance 
with respect to a — > a' = ae dl , response functions for the </>± -field are derived as 
(Appendix A) 



^sin^(r) = 2(T T sin^_(x,r)sin0_(O,O)) 

r ln(r/a) l 



exp 



di 



• exp 



+ 



Rcos 4>- (r) = 2 (T T cos 0_ (x, r) cos (0, 0)) 

r\n(r/a) ^ 



exp 



Jo 



di 



■ exp 



ln(r/«) 



r\n(r/a) 
J o 92cf>-(l)dl 



g 2( f > -(l)dl 



(3-5) 



(3-6) 



R S m<t> + (r) = 2 (t t sin ^ + (x,r)-^-x^ sin0 + (O,O)^ ^/ cos ^-^-x 



exp 



ln(r/a) 



• exp 



ln(r/a) 

- / g2<t>+{l)dl 



f 

Jo 



(3-7) 



^cos 0+ (r) = 2 ( T r cos ^+ (x, r) 

^ ln(r/a) 



2t \ , A / (2t 

— x cos 0+(O, 0) ) / cos — . 

VF J I \VY 



-x 



exp 



Jo 



K<t,(l)dl 



■ exp 



rln(r/a) 

/ 924>+{l)dl 
Jo 



(3-8) 



where the term with the quantity tan _1 (vFT/x) has been neglected for the present 
numerical calculation with small t. In Eqs. (3-5)^(3-8), the quantities K,p(l), g 2( j > - (0 
and g2(fi+ (0 are calculated from the renormalization equations, 



= -\ 9 % + {l) Klil) Jo (-^ ) + tgl,_(l) 



dl 



g 2<p+ (l) = (2-2K (j> (l)) g 2<t>+ (l), 



(3-9) 
(3-10) 
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J l 92< t >-{1) = (2 -VKS)) 92<p-(l), (3-H) 

±(^)=^-^WKMJ 1 (<M), (3,2) 

where J u (z) is the v-th Bessel function. Equations (3-5)~(3T2) have been derived by 
use of a method similar to that of Giamarchi and Schulz. The initial conditions 

inEqs. (3-9)~(3-12) are chosen as iQ,(0) = 1, 52^+ (°) = "(52-52)) 520- (0) =92-92 
and £(0) = i. We note that the renormalization equation (3- 12) for small t becomes 
equal to that derived from the perturbation in terms of t.ffl The negative sign of the 
second term of Eq. (3- 12) 03 indicates the reduction of interchain hopping. 
From Eqs. (|3-2|), (3-5)~(3-8), the response function can be expressed as 



R A (r) = R a (x,t)R a {x,t), (3-13) 

where the largest response function and the next one in R A (r) with fixed length 
r(= \Jx 2 + (vft) 2 ) (i.e., the corresponding energy vp/r) lead to the phase diagram 
for the dominant state and subdominant state, respectively. 

In the present paper, we assume that the response function R A (r) is equivalent 
to that of the Fourier transform with u = vp/r. The validity of such a treatment is 
discussed in S5. 



§4. Response functions and phase diagram 



We examine the renormalization 
group equations given by Eqs. (3-9)~ 
(3T2), where parameters are chosen as 
t/ivpa- 1 ) = 0.0025 and g 2 - g' 2 = 0.02. 
The /-dependence of K^,(l), g2 < j>-(l) and 
924>+{l) are shown in Fig. 2. Increas- 
ing I, K^il) and g 2( f > -(l) increases to the 
strong coupling regime, while #20+ (0 re- 
duces to zero. Note that the quanti- 
ties d(l) of Eq. (A-8), which is ob- 
tained by substituting the results of 
Eqs. (3-9)~(3T2), is invisible in the 
scale of Fig. 2. The actual vari- 
ation of g2<t>- (I) and 520+ (0 appears 
for l> ln[vFa -1 /4t] ~ 4.6, which is 
designated by the arrow. This result 
indicates the validity of the previous 
treatment IIP that and <?2<£± have 
been replaced by the initial values for 
e l < vpa~ l /t. Since the Bessel func- 
tion Jo(4t(/)/(wFa -1 )) becomes small 
for e l > VFa~~ 1 /t, the 0+-term can be 




|5 

I 



Fig. 2. Quantities K^ (l) , g 2< t,- (I), £20+ (0 and 
At{l) / {v-poT 1 ) as functions of I — ln(r/a) 
for t/^pa" 1 ) = 0.0025 and g 2 - g' 2 = 
0.20. The arrow denotes I = lo = 
ln[u F a _1 /4i] (- 4 -6)- 
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neglected and leads to strong coupling 
due to the 0_-term. Equation (3-12) 
reveals the fact that the hopping en- 
ergy is reduced by the mutual interac- 
tion. Therefore the range for the one- 
dimensional regime is examined in de- 
tail. 

In Fig. 3 the present results (solid 
curves) are compared with the previ- 
ous ones (dotted curves) , which were ob- 
tained by using initial conditions, such 
that the terms including t(l) are dis- 
carded in Eqs. (3-9) and (3-12) and 
for which K(Iq) = 1 and g2<p- r lla) = 
92 ~ 9' 2 with l = ln[u F a- 1 / 4 *]-i)'IZI ) In 
the solid curve of K^(l) — 1, the tan- 
gent, which displays the linear depen- 
dence for 1<^Z<^5, exhibits a rapid vari- 
ation for I ~ 6, indicating the crossover 
of the scaling property from the one- 
dimensional regime to the regime of rele- 
vant hopping. Therefore we define such 
a point, I = li, by the condition that 
the absolute value of the variation of 
the tangent becomes maximum. Good 
agreement between the solid curve and 
the dashed curve is obtained for I > h, 
where l\ > hi^pa /4i]. For the small 
I (< h), Eqs. (3-9)~(3-ll) up to the 
calculated as (Appendix B) 



10' 



10" 



1 1 1 1 


i i i I 


/ i i i i : 


/ (82-82) 
1 1 1 1 1 1 



10 



/ 



Fig. 3. The Z-dependence of K${1) — 1 and 
g24>-(l) — (<?2 — 172)1 which denote the vari- 
ation from the initial values. The param- 
eters are the same as in Fig. 2. The dot- 
ted curves represent the previous results, uP 
which were obtained by use of the initial 
condition that Q2<t>+ {k)) = 0, K(lo) = 1 and 
ff2</>- (lo) = 92 ~ g'2 with l = ln^Fa" 1 /4t]. 
The factor 2 error in the renoHnalization 
equation of the previous papers™'™ is cor- 
rected. The quantity fa denotes I for the 
onset for the relevant behavior of K<f,(l) in 
the presence of t. 



lowest order of t/(v-pa x ) and #2 — g'2 are 



' ^2 



K4.ii) = 1 + (52 - g' 2 ) 



VpCi 



e 2l -l) + 



WO = -(52 - g'2) + - g' 2 f f-rzr) ( 

92<j>-{l) = (52 - g' 2 ) + (52 - g'2) 



,21 



1-21) + 



J \3 ( * 



1-2Z ) + ■ 



(4-1) 
(4-2) 
(4-3) 



Equations (4-1) and (4-3) can reproduce the numerical results sufficiently, i.e., 
the linear dependence of the solid curves for I < l\ in Fig. 3. Thus it turns out that 
there is a perturbational effect of t for the one-dimensional regime which is obtained 
for I < l\. The dependence of the crossover energy on the interaction is discussed 
later. 
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By substituting the solutions of Eqs. 
(3-9)~(3-12) into Eqs. (3-5)~(3-8), we 
obtain response functions for <f)± which 
are shown by the solid curve in Fig. 4 
for t/(v F a~ l ) = 0.0025 and g 2 - g' 2 = 
0.2. The dotted line corresponds to 
those in the absence of the hopping. 
For ln(r/a) < Zi(~ 6), the response 
functions of R sin( j } _(r) and i? C os </»+(/) 
(-Rcos <£_ ( r ) and ^^^(r)) are almost 
the same as that for t = (dashed line), 
and their deviations from the dashed 
line due to the interchain hopping be- 
come visible for ln(r/a) > l\. Since the 
renormalization group method is based 
on perturbation with respect to g2<f>+ 
and g24>- , the present result is valid only 
for the regime of weak coupling. For 
large ln(r/a), -R S in</>_(0 increases due 
to the strong coupling but is actually 
reduced to a finite value. We define 
ls(= ln(r3/a)) as the location for the 
crossover from the weak coupling regime 
to the strong coupling regime which cor- 
responds to a minimum of R S m<f>^( r )- 
As is shown in Appendix C, the exci- 
tation gap, A, in Eq. (2-4) can be es- 



timated from A = v-pa 
2(K 4> (l A ) - 1) ~ 1, i.e., 



exp[— Ia\ with 



A = At exp 



10' 



10 



10 



10 



-6 



~i 1 1 1 1 1 1 r 



sin(# + 



5 10 

!/(&-&') 

_l I I I I I I I I I I L_ 



5 10 

ln( r I a ) 

Fig. 4. The response functions of the <f>±- 
field as a function of \n(r/a), where 
r = \/ x 2 + (vft) 2 . The quanti- 
ties R s in<p_{r), R C os<p_(r), R s in<p + (r) and 
-Rcos 4> + (r) are defined by Eqs. (3-5), (3-6), 
(3-7) and (3-8), respectively. The parame- 
ters are the same as in Fig. 2. The dotted 
lines denote the response functions with 
t = 0, and the location, h(= ln(r3/a)), 
corresponds to a minimum of -R S m,£_M- 
In the inset, ln[7? s i n </>_ (t^)] (solid curve) as 
a function of g2 — g'2 is compared with the 
analytical one (the dash-dotted curve). 



7T 



9*2 



+ 1 



(4-4) 



In the inset in Fig. 4, lni? s i n ^_(r3) (solid curve) as a function of 1/(52 — g'2) 1S 
compared with the analytical function -R s in0_( r zl) (dash-dotted curve) where I a = 
hi^zi/aO- Since ~ r/\, it is found that the minimum of i? S m0_ ( r ) at ^3 corresponds 
to the formation of the gap. 

Now we examine the actual response function of Eq. (3T3), which is a product 
of the response function of the total fluctuation and that of the transverse one. 
In Fig. 5, the response functions with the largest (the dominant one), the second 
largest (the subdominant one) and the third largest in magnitudes are shown for 
t/vpa' 1 = 0.0025, g 2 - g' 2 = 0.20 and g 2 + g 2 = 0.10. In previous paper, the 
dominant and the subdominant states were examined in the limits of both long range 
(Fig. 1(b)) and t — * 0. In the present study, by tracking a state with arbitrary length 
scale (corresponding to the inverse of the energy) we found that the crossover of the 
second dominant states from Cm ; n to S_im occurs at m(r/a) = l 2 . In the region of 
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^ 10 
let; 



10 



f/( VpCt ) = 0.0025 
g 2 -fo' = 0.20 
g 2 + g 2 - = 0.W 



S" ] ,in~~ *" 



_l I I I L. 



_l L 



5 10 

ln(r/a) 

Fig. 5. The normalized response functions, 
Ra{t), of the order parameters, C|| jOUt , 
C|| jin and yS±,in as functions of ln(r/a), 
where t/(u F a _1 ) = 0.0025, 92-92= 0.20, 
32 + 32 = 0.10 and notation is the same 
as in Fig. 1. The dotted line denotes the 
response function in the absence of inter- 
actions. With increasing ln(r/a), the dif- 
ference between C|| j0ut and C|| >in becomes 
notable for ln(r/a) > h, and a subdomi- 
nant response function varies from C|| jin to 
Sj_,i n at ln(r/a) = h. 



15 



r/OpOC ) = 0.0025 
g 2 -g 2 ' = 0.20 




ln( v F a / t ) 



0.1 

82 + 82' 



0.2 



Fig. 6. Quantities h and h as functions of 
92 + g'2, with fixed t/^pa" 1 ) = 0.0025 and 
92 — g'2 = 0.20, where the one-dimensional 
regime, i.e., the irrelevant t(l) is obtained 
for < ln(r/a) < h and the subdom- 
inant state is given by Ciu n (Sa,i n ) for 
\n(r/a) < h Qn(r/a) > h). The value l 3 
(dashed line) is the upper bound for the 
present treatment of the renormalization 
group method. The maximum value for 
the one-dimensional regime denoted by h 
is larger than \a\vYa~ 1 /At] (dotted line). 



ln(r/a) < li, there is still a small effect of the interchain hopping which removes the 
degeneracy of C_i_ iOUt and C|| ; i n . Such an effect can be seen analytically: Eq. (3-13) 
with m(r/a) <C h is expanded as 



^C||.out( r ) 



a\ Ke fa 



1 + 



l-(92-3 2 ) 



r / \r 

n - g' 2 ? ( 4t(0) ^ 2 



32 



e 2Mr/a)_ 1 _2in[- + 







a\ Ke fa\ 1 ~(f2-9 2 ) 



,M2 



(7) 



q, N l-(ff2-fl 2 ) 

r 



1 + 



(fl2-^ 2 ) 2 /4*(0) 



32 



-1 



e 21n(r/a) _ 1 _ 21n 



n 



+ 



, (4-5) 



, (4-6) 



(4-7) 
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For the range satisfying ln(r/a) > I3 corresponding to an energy lower than A, there 
appears a gap in the transverse fluctuation due to the relevance of the interchain 
hopping which leads to a large enhancement for Cm out and iSj_ )0 ut an d strong sup- 
pression for C11 j n . We have found two kinds of quantities, l\ and I2, where coupling 
constants are renormalized toward the relevant t in the region ln(r/a) > h and S± i i n 
is found as the subdominant state due to the relevant t in the region ln(r/a) > h- 
In Fig. 6, l\ and I2 are shown as a function of 52 + 92- The dashed line denotes £3 
corresponding to the upper bound of ln(r/a) for the present treatment of the renor- 
malization group method and the dotted line denotes m^pa -1 /*]- In the region 
with ln(r/a) < l\, the state exhibits one-dimensional properties since the interchain 
hopping can be treated perturbatively. We note that I2 < h for small §2 + g'2 due 
to I2 ~ (§2 + 92){92 — 92)~ 2 {vfoi~ 1 /t) 2 which is obtained from Eqs. (4-5), (4-7) and 

{92 +92)1 < 1- 

Based on the results of Figs. 5 
and 6, we obtain the phase diagram as 
a function of 52 and g' 2 for the domi- 
nant and subdominant states with fixed 
ln(r/a). The boundary between Cm in 
and S± t i n is obtained from Fig. 6 
by estimating 52 and g' 2 which satisfy 
ln(r/a) = I2 with fixed r. For exam- 
ple, the subdominant state, S± j j n , with 
ln(r/a) = 12 appears in case 92 — 92 = 
0.20 and g 2 + g' 2 < 0.11 due to the 
enhancement of the charge fluctuation. 
We note that the response function with 
fixed r corresponds to the response func- 
tion for finite energy with uj = v-p/r. In 
Fig. 7, the phase diagram in the 52-^2 
plane is shown with t/(vYOt~ l ) = 0.0025 
and fixed ln(r/a) = 12. The solid line 
denotes the boundary for the dominant 
state and the dotted line is the bound- 
ary for the subdominant state. The re- 
sults of Figs. 5 and 6 correspond to the 
state in the region #2 > Ifl^h i- e -> and 
Ih, which show the states with low en- 
ergy and high energy, respectively. The 
phase diagram of the entire region has 
been obtained by use of the symmetry property of order parameters which is given 
by Eqs. (2-8)~(2-ll). The region of J L (J = I, II, III and IV) which is similar to 
the state in the zero limit of energy decreases rapidly by the decrease of ln(r/a). 
We have found that the region of Jl becomes very narrow in case ln(r/a) ~ l\ 
(~ ln^poT 1 /*])- Therefore the phase diagram is similar to Fig. 1(a) for most of the 
region in the <?2-<72 plane in the case that ln(r/a) < l±, i.e., for energy larger than t. 
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Fig. 7. Phase diagram for several kinds of or- 
der parameters in the plane of <?2 and g' 2 , 
where the dominant states (subdominant 
states) are given by C||, out (S±,i n ), C|| >out 

(C||,i n ), C±,out (Sll.in), C"±,out (C±,in), S\\j n 

(C±,out), •S'n.in (<5|| >0 ut) and S±j n (Cy, out), 
S±,in (S'x.out) for the regions II, Ih, IIl, 
IIh, III h and IVl, IVh, respectively. 

The dotted curve, which corresponds to the 
boundary for the subdominant states is cal- 
culated for ln(r/a) = ln(u F a~ 1 /u;) = 12. 
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Finally we examine the characteris- 
tic energy u\ which separates the one- 
dimensional regime from the relevant 
regime of t. The length corresponding 
to u>i is shown by l\ = ln^pa" 1 /lo\] in 
Fig. 3. In Fig. 8 the quantity uj\jt is 
shown function of t/(vpa 1 ) with 
fixed g 2 -g' 2 = 0.1, 0.2, 0.3 and 0.4. The 
quantity to\/t decreases by the increase 
of 92 — 92 an d increases with increasing 
t/ {vpa~ l ). The dotted line which is ob- 
tained by use of the extrapolation of the 
solid curve is expressed as uj\ jt = const • 
[t/ (fFa _1 )] ai , where a\ depends on g 2 — 
g' 2 . Such a characteristic energy lo\ can 
be estimated from Eq. (3-12) by as- 
suming that uj\ ~ u[ = vpa -1 exp[— l[] 
with til'^/ivFa- 1 ) = O(l), where l[ ~ 
l\. Since Eq. (3-12) is calculated as 
£(7)/0Fa -1 ) = t/ (vFa" 1 ) exp[{l - (g 2 - 
g' 2 ) 2 /2}l] for t/(v 7 a' 1 ) < 1, we obtain 




f/(v F oc ) 

Fig. 8. The ratio of uji to t as a function of 
t/(upQ _1 ), where lui = upa -1 exp[— h]. 
The dotted line is the extrapolation which 
is given by using the equation ui/t oc 
(t/tiFd _1 ) ai . In the inset, ot\ is shown as a 
function of (72 — g'2, where the dash-dotted 
curve denotes Qi = (\/2)(g2 — g'2) 2 ■ 



const • t 



(t/ivpa- 1 } 



(92-g' 2 ) 2 /2 



(4-8) 



In the inset of Fig. 8, ct\ obtained from the data of the solid curve is represented 
by the closed circle which is compared with the a\ = (g 2 — g 2 ) 2 /2 of Eq. (|4-8| ) 
(dash-dotted curve) . The good agreement of the two results indicates the validity of 
Eq. ( f4-8| ) as a characteristic energy u\. Here we note the assertion of BourbonnaisEp 
that one-dimensional fluctuation does exist at temperatures lower than the energy 
of the transverse hopping t± in quasi-one-dimensional systems. The suppression of 
the crossover temperature T x (< t±) in a quasi-one-dimensional system comes from 
the reduction of the density of states at the Fermi energy, which is caused by the 
mutual interaction. u> We find that u\ of the present calculation is essentially the 
same as T x since a system of two-coupled chains already shares common features 
with a quasi-one-dimensional system regarding the role of interchain hopping. 



§5. Summary and discussion 



In the present paper, we have examined the role of interchain hopping in the 
entire energy region in two-coupled chains of the Tomonaga model by use of the 
renormalization group method. From the calculation of response functions for CDW 
and SC states in the real space, we obtained the phase diagram in the plane of 
mutual interactions in which the dominant and subdominant states are shown with 
fixed energy or temperature. The interchain hopping exhibits a perturbational effect 
for energy higher than u\ (= v-poT 1 e - ' 1 ), while the degeneracy of in-phase and out- 
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of-phase parings is removed due to the relevant behavior of the hopping for energy 
lower than u\. The hopping leads to the formation of the excitation gap in the 
transverse fluctuation of the 0-field for energy lower than vpa -1 e~ lz . The crossover 
energy, uj\, is lower than t due to the interaction which gives rise to the reduction of 
the density of states at the Fermi surface. 

We have examined the two-chain system by use of the renormalization group 
method in real space in order to examine the states at finite energy or at finite tem- 
perature. We comment on the relation between the present result and that derived 
by the scaling of the energy cutoff, i.e., in Fourier space. FabrizioliS has studied 
the same model by applying the renormalization group method with the scaling of 
energy to Eq. (2-1), where all interactions are calculated perturbatively. In the 
present calculation, the two-chain system is transformed into the phase Hamilto- 
nian through the bosonization method, where only the non-linear terms #20- and 
g2<p+ are treated perturbatively. By replacing l{= ln(c//a)) with — ln(u' /u}o), where 
wo(~ v-f/oc) is the band cutoff, it is found that the renormalization equations of the 
coupling constants obtained by Fabrizio are the same as ours, Eqs. (3-9)~(3T2), up 
to the lowest order of the coupling constants and the interchain hopping. But there 
is a difference in the terms of higher order in the interchain hopping. Invariance of 
the renormalization equations with respect to t —* — t exists in the present result, 
Eq. (3T2), but does not for his results. It is not yet clear if such a discrepancy 
comes from the present choice of the phase Hamiltonian or from some other source. 

Finally, we discuss the response functions at finite energy. Since, to our knowl- 
edge, there is no calculation for response functions of two-coupled chains in terms 
of the renormalization group method with an energy cutoff, we mention the case 
of no misfit parameter, i.e., a one-dimensional Fermi gas. The response function 
for the energy cutoff is calculated by applying the renormalization group method 
to the quantity given by Ri(0,u>) = ttvf (d Ri(0, u>)/d In uj),i^ where Ri(k, u>) is the 
Fourier trans for m of Ri(r). By comparing Eqs. (3-5)~(3-8) in the case t = OliHP 
with Ri(0,u>),t3 it turns out that Ri(r)/(a/r) 2 with r/a = ojq/uo can be replaced by 
Ri(0,uj) when the /-dependences of K^il) and g2<j>±ij) are neglected. In the present 
calculation, such a replacement is justified for Z<10 and is meaningful qualitatively 
for larger /. Thus, the present results may correspond to those at finite energy or at 
finite temperature. 
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Appendix A 

Derivation of Renormalization Equations 



The response functions are calculated by use of the renormalization group method 
as follows. HIP Treating the nonlinear terms in Eq. (2-4), i.e., g2(j>+- and g20--terms, 
as the perturbation, the response function for R s i n j,-(x, r) of Eq. (3-5) up to the 
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second order is calculated as 

R S m<t>-{xi -x 2 ,n -r 2 ) = 2(T T sin<^_(£i,Ti)sin0_(z 2 ,T 2 )) 



exp 



--r^U(n - r 2 ) 



dr /r\ 3 ~ 4K - 



d 2 r 3 



e -(2/^)tf(n-r3) e -(2/J^)tf(ra-rs) 



Cx 



2 

92<f>-\ 



2 

520^ 



a \a 



At 



Jo [ —r ) U(n - r 2 ) 
•^2 ( — r I cos2# ri _ r „ 



1 o 1 f dr / r \ 3 " 4 /^ 



+ ^^20-^2 / — l~ 



l/(ri - r 2 ) + 



(A-l) 



where r = yjx 2 + (vft)' 2 , £7(r) = \n(r/a) and tan# r = v-pr/x. J v [z) is the z/-th 
Bessel function. By raising the expansion of r.h.s. of Eq. (A-l) to the exponential 
and using a as the lower cutoff of the integrations, one obtains 



R^<i>-{xx - x 2 ,n - r 2 ) 
l 



exp 



U(n - r 2 ) - d cS cos(20 ri _ r2 ) 



1 + —g u _ e (2/K^)U( ri -r 2 ) f e -(2/J^)tf(ri-r 3 ) e -(2/K^)U(r 2 -r 3 ) + . 

47T * 7a a 2 



(A-2) 



where the effective quantities are given by 



1 



dr f r 



3-4fO 



1 2 



2^-J a a \a) 
■h —r I . 



, (A-3) 
(A-4) 



The parameter d is introduced in the zeroth order, where d = at I = 0. By assuming 
the scaling relations that Eqs. (A-3) and (A-4) are invariant for a — > a' = ae*, the 
renormalization equations for the coupling constants are obtained as 

|^(0 = -^.(0*2(0 Jb (^r) + ^ 2 V(0> (A-5) 



^ + (/) = (2 -2^(/))^+(0, 



(A-6) 
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-g 2 ^(l) = (2-2/K 4> (l)) g 2 ^(l), 



dl 



d(l) = -9^(1)^(1) J, 



( m 



(A-7) 
(A-8) 



where a(0) = a in the Bessel functions. The renormalization for At{l) / (vpa -1 ) is 
calculated as follows. The actual difference of the density between the two bands of 
Eq. (|2-2j), An, is given by 



T f 

An = —2ta/vp + — a I dxdr (d x 4>+) , 



(A-9) 



where d x (f> + = T^S p (^r ±4> p ,+ — ip p _tp p —). The renormalization equation for t(l) is 
derived by assuming the scaling relation for Eq. ( |A-9j ). The second term on the r.h.s 
of Eq. ( |A-9j ) is rewritten as 



dx dr(d x ( 1 



'+/ 



1 

Z 



Tr 



8 

ox exp 



dx dr d x (j> + e- S dTH 
-J drTC + X J dx dr d x <p + 



A=0 



(A-10) 



where Z = Trexp[— J dr7i\ . By writing 
Eq. (A-9) is calculated as 

2t 



<j) + - (2Tr/v F )K <t) Xx in Eq. (A-10), 



An 



VpCt 

2t 
vpa~ 



2 T f I ( At 

i + -92<t>+K (j) — J dxdr (xsm. f 20 + - — x 



oo ^ r 

, a 



2-4K 



■h ( -r ) + 

Vp 



The infinitesimal transformation a' = e dl a for Eq. (A-11) leads to 



d_ 

dl 



-i 



VpOt 



vpa 



g% + {l)K^l)J, 



\VFa~ 1 



(A-11) 



(A-12) 



Note that the present model is similar_to but slightly different from the model of a 
single chain of fermions with spin 1/2. IliP The hopping t corresponds to the external 
magnetic field in the latter model, where the renormalization equation for the mag- 
netization has been derived by assuming that the external magnetic field satisfies 
the scaling relations. 

The response functions are calculated by writing Ri(x, r) = Fi{r) ■ exp[— K^ 1 In 
(r/a) — dcos28 r ], where i = sin <p_ , cos 4>- > sin <p + and cos<p + and the ± sign cor- 
responds the response function for the <p±-&eld. By assuming the scaling relation 
F i {r,a(l),K <p (l),g 2 $±(l)) = h (dl, K^(l), g 2 <p±(l)) ■ Fi(r,a(l + dl), K^l + dl), g 2(t>± (l 
+dl)), the multiplicative factor l{ for i^ s i n <A- is obtained as 



Ism<f>-{dl, Ktf,, g2<j>±) 
1 



exp 



924>-dl - -92<t>+ J v 



At 



vpa 



-i 



U(n - r 2 )dl 
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At 



4%>+ J 2 ( - — x j cos 26 ri - r2 dl + -g^j^Uin - r 2 )dl 



(A-13) 



Thus the functions F{ are expressed as© 
Fi(r, Kf, g 2 <f,±) = exp 



ln(r/o) 

£ Hl i {dl,K 4> {l),g 24>± {l))] 
1=0 



(A-14) 



Equations ( |A7T2| ) (A- 13) and (A-14) lead to 



exp 



ln(r/o) 



dl 



1 



^(0 



cos 26 r 



(A-15) 

where other response functions are obtained in a similar way. We note that the last 
term of Eq. ( [A- 15 ) is discarded in Eqs. (3-5)~(3-8), i.e., Ra(t), due to the condition 
d <C 1 in the present numerical calculation. 

Appendix B 

Derivation of Eqs. (4'1) ~ (4'&) 



We examine the renormalization equations for small Noting that 4t(l)/(vF 
a -1 ) <C 1 for the small l/h, we use the approximation that Jo(4t(/)/(wpa -1 )) ~ 
1 - {4t(l)/(v F a- 1 )} 2 /4 and Ji(4t(Q / '(u F a -1 )) ^ {^/(vpa -1 )} /2. By retaining 
terms up to the order of (g 2 — g' 2 ) 2 , Eq. (3-12) is calculated as 



t exp 



(Bl) 



where t = 4t/(vpa *) and g = g 2 — g' 2 . By rewriting K^l) as K${1) = 1 + go(l)/2 
and making use of Eq. (B_l), Eqs. (3-9), (3-10) and (3-11) are calculated as 



/ 2 



flb(0 = ^ 2 (e 2 ' - l) , 
924,+il) = -9 + {^9 3 (e 21 -1-21' 
520 _(O = <7 + ^ 3 (e 2 '-l-2Z), 
which correspond to Eqs. (4-1), (4-2) and (4-3), respectively. 

Appendix C 

Estimation of the Energy Gap 



(B-2) 
(B-3) 
(B-4) 



We examine Eqs. (3-9)~(3T2) with small coupling constants for I 3> h, where 
the term including J v (4t(l) / (vpa" 1 )) and <?2</>+ can be disregarded. In this case, the 
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renormalization group equations are rewritten as 



(C-l) 
(C-2) 



where the initial condition is given by K^Iq) = 1 and g2<f>-{h) =92 — 92 with 
Zo = — ln[4t/(fFO -1 )]- The differential equations (C-l) and (C-2) can be integrated 
as 

<&-(0 - 8 W) - k#*(0) = ^-(io) - 8(^(i ) - Ib^(Io)) 

* 5 2 V(0 " flg(0 - 8 = (52 - 3 2 ) 2 - 8, (C3) 

where K^(l) = 1 + go(l)/2. From Eqs. (C4)~(C-3), one obtains 

9o(l) = \92 ~ 92\^[\92 ~ g' 2 \(l ~ k)] • (C-4) 
By defining the gap A as A = upa -1 exp[— Z^i] with </o(^d) = 1> on e obtains for A, 

IT 1 



A = 4t exp 



2 |52 - 52 



+ 1 



(C-5) 



From Eqs. (C-3)~ flC-5|) the response function -Rsin (r) at Ia{= hi(r^\/a)) is calcu- 
lated as 



(r A ) = exp 



dl 



1 



X ( 1 + 



o V 

4^ \i-tea-$£) 

exp 

1/4 



7T 1 



(32 "3 2 ) 2 



2 1 32 - 92 
1 



7T + 1 



132 - 3 2 



77 + J1 + 



(32 - 3 2 ) 2 



,(C-6) 



which is shown as a function of 1/(32 — 3 2 ) in the inset in Fig. 4. 
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